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The Demise of the Classical BLR in the Luminous Quasar PG1416— 129 
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ABSTRACT 



New observations of the broad-line quasar PG1416-129 reveal a large decline in its 



o 
o . 

^^ , continuum luminosity over the past ten years. In response to the continuum change the 

"classical" broad component of H/3 has almost completely disappeared (a x 10 decrease 



in flux) . In its place there remains a redshifted/redward asymmetric very broad emission 
line component. The significance of this change is multifold. (1) It confirms the existence 

^^ ■ of a distinct redshifted Very Broad Line Region (VBLR) component that persists after 

the demise of the broad component and that is frequently observed, along with the broad 

,__i , component, in radio-loud sources. (2) The smaller (x2) intensity change in the H/3 very 

broad component supports the previously advanced idea that the VBLR is physically 

£\j | distinct and likely to arise in an optically thin region close to the central source. (3) 

The presence of a strong very broad component in the radio-quiet quasar PG1416- 
129 reinforces the notion that such "population B" quasars share similar spectroscopic 
(and hence geometrical and kinematical) properties to radio-loud sources. (4) AGN 

-^ . can show broad, very broad, or both line components simultaneously, making statistical 

comparisons of source profile widths difficult. (5) The interpretation, in reverberation 
JL \ studies, of the presence or lack of correlated response in broad line wings will be affected 

H \ by this composite BLR/VBLR structure. 

cd - 

Subject headings: quasars: emission lines -- quasars: general -- quasars: individual 
(PG 1416-129) — line: formation — line: profiles 



1. Introduction 

We present new spectroscopic observations of the radio quiet quasar PG 1416-129 that reveal 
a significant decrease in optical continuum intensity and the virtual disappearance of the "classical" 
broad component of H/3 (hereafter indicated as H/?bc) last observed 10 years ago. There remains 
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a strong, very-broad (indicated as H/3vbc) an d redshifted component that has declined much less 
dramatically in response to the continuum change. In §2 we describe the new observations, data 
reduction procedures and a comparison of the old and new spectra for this source. This is followed 
(§3) by a discussion of physical and phenomenological implications of the spectra change in PG 
1416-129. 



2. Observations and Analysis 

New spectra of PG 1416-129 (z=0.12928±0.00007 from [OIII] A5007) were obtained on the 
2.12-m telescope of the Osservatorio Astronomico Nacional at San Pedro Martir, Baja California, 
Mexico on May 9 2000 between 7:28 and 8:32 UT, under photometric sky conditions. A Boiler &: 
Chivens spectrograph was equipped with a 600 1/mm grating with a Tektronix TK1024AB CCD 
chip as the detector (2 arcsec slit width). Suitable standard stars were observed at a similar 
airmass. The spectrum covers the range AA4438 — 6537A with an effective spectral resolution of 
3.5 A FWHM. It was reduced with standard IRAF procedures. Line properties were derived after 
normalization of the two spectra to the [Oiii] AA4959,5007 flux in a spectrum obtained on February 
17, 1990 (Boroson & Green 1992: hereafter BG92). 

Figure 1 shows the 4200-5500A rest frame region in both the 1990 and 2000 spectra. After 
subtraction of the underlying continuum (fit with a low order polynomial; contamination from the 
light of the host galaxy is negligible), the H/3 profiles were cleaned of [Oiii]AA4959,5007, H/3nc> 
and HeiiA4686 emission (for details see Marziani et al. 1996). Observed line flux, rest-frame 
equivalent width and FWHM values derived from the 2000 May 8 spectrum are given in Table 
1. Errors are 20% (2a) for the flux and equivalent width and 10% (2a) for the FWHM. FEll opt 
emission is weak in the 1990 spectrum (BG92 estimated W(Feii A4570)= 31 A) and weaker in the 
new spectrum as expected if it correlates with the strength of the classical BLR (e.g. BG92). We 
subtracted a scaled and broadened I Zw 1 template (-Fa <1x 10 -14 ergs s -1 cm -2 ) from the 2000 
May 8 spectrum in order to remove FEii opt contamination. HeiiA4686 emission may have a more 
significant effect on the blue side of the H/3 profile if a significant HeiiA4686 VBC exists in our 
new spectrum. After several empirical attempts, the smoothest residuals were obtained by using 
the H/3vbc measured from our new spectrum as a template to subtract suspected broad HeiiA4686 
emission. This subtraction assumes a much stronger HeiiA4686 emission than the estimate given 
in BG92 (W(HeiiA4686)/W(H/3) ss 0.04). The result of a spline fit to H/3 V bc after HeiiA4686 
subtraction is shown for both spectra in Fig. 2 with parameters detailed in Table 2. 

The main observational results of this investigation are: (a) the almost complete disappearance 
of the H/3bc ( on 2000 May 8 the H/3bc A ux is ~ jn OI " the original flux) in a reasonably luminous 
quasar with bolometric luminosity ~10 45 ergs s _1 and (b) the much smaller variation in the H/3vbc 
component (a factor ~ 2). 



3. Implications 

3.1. Implication 1: Existence of a VBLR Associated With H/3vbc 

The 1990 observation of the H/3 region in PG 1416-129 showed inflections between narrow 
(H/?nc)j broad and very broad line components. The 2000 spectrum shows the almost complete 
disappearance of the broad component H/3bc made more evident by the lack of change in H/3nc- 
The bulk of the residual broad line emission is associated with the very broad component H/?vbc > 
also seen in 1990. The robustness of the H/3vbc component to a dramatic continuum change that 
quenched H/3bc strongly suggests that the VBC represents a distinct emitting component. Similar 
arguments have been made previously based on variability and line inflection evidence (Corbin 
1997a,b). The weakness of FEii opt emission in PG 1416-129 suggests that the H/3vbc cannot be 
attributed to a red-shelf of enhanced FEH op t emission (Korista 1992). Do all RL and RQ (pop B: 
with FWHM(H/3) > 4000 km s^ 1 see Sulentic et al. 2000a) sources show a VBC? The answer is 
"no", at least for a VBC like the one observed in PG 1416-129, within current detection limits. 
Sources exist with no hint of a BC/VBC inflection, H/3 red wing, or profile shift. 

Evidence also exists for VBC components of HeiiA4686 (e.g., Ferland Korista & Peterson 1990; 
Marziani & Sulentic 1993) and CIVA1549 A (Laor et al. 1994; 1995) in several sources. It is very 
likely that these VBC components are related to the H/3vbc> although HeiiA4686 may often show 
a more symmetric, boxy profile (as in PG 1138+222; Marziani & Sulentic 1993). However, even if 
it is possible to disentangle the HeiiA4686 emission from the surrounding FEii opt contamination, it 
is usually weak, with a flux distributed over a wide wavelength range often merging smoothly with 
the blue wing of H/3. HeiiA4686 asymmetry measurements are therefore intrinsically difficult and 
dependent on assumptions about the H/3bc shape even if the s/n of the data is high. Note that 
the classical BLR of CIVA1549 and other high ionization lines is called VBLR in the ILR/VBLR 
scheme followed by Brotherton et al. (1994) (see Sulentic & Marziani 1999). In Corbin (1997ab) 
ILR= our NLR+BLR but VBLR definition is the same. The CivA1549 VBC that we refer to is 
redshifted and much broader than the BC in analogy to H/3 (see e.g. Laor et al. 1994, 1995). 



3.2. Implication 2: Physics of the BLR and VBLR 



1 12 



BLR The light travel time tlt for a spherically symmetric BLR is tlt ~ 500L^ 42 f f j^nj^ days, 
where ff,-7 is the filling factors of the line emitting gas in units of 10~ 7 , Ljj/3,42 is the H/3 luminosity 
in units of 10 42 ergs s -1 , and n is the number density in units of 10 11 cm~ 3 . The H/3bc luminosity 
of PG 1416-129 is L h/3bc w8x 10 41 ergs s" 1 (excluding VBC emission; H = 75 km s _1 Mpc _1 , 
qo = 0.5). The tlt for the BLR of PG 1416-129 is therefore most likely between several months 
and a few years. The almost complete disappearance of the H/3bc i n ~ 10 yr does not pose any 
fundamental challenge. The simplest scenario is that H/3bc faded following a high amplitude (a 
factor ^ 4) decrease in ionizing continuum level. The H/?bc luminosity can be conventionally 



explained as due to an ensemble of optically thick clouds of very small filling factor (ff ~ 10 7 ) in 
a spherically symmetric region extending from the VBLR radius up to 10 19 cm. 

Another potentially very important constraint arises from the apparent absence of a BC in the 
HeiiA4686 profile of PG 1411-129. Examination of the PG spectral atlas (BG92) indicates that 
when a broad component of HeiiA4686 is visible, it is usually broader than H/3bc- This is especially 
true for population B and RL sources (e.g. Corbin & Smith 2000). In PG 1416-129 there is evidence 
for a HeiiA4686vbc but only in our new spectrum where H/3bc has almost disappeared. Negligible 
HeiiA4686 emission (I(HeiiA4686bc)/I(H/3bc) ~ 0.03) implies that the ionization parameter must 
be very low, T ^ 10 -4 in the BLR of PG1416-129. This condition is satisfied if the emitting gas 
is either high density (n c ^ 10 12 cm~ 3 ) or located farther away from the continuum source (a 
condition not supported by the line width and by the variability timescale). If most HeiiA4686 is 
produced in a VBLR rather than a BLR, the problem of an ionizing photon deficit (Korista et al. 
1997) may vanish and we may need to reconsider the formation of the high ionization lines (HIL) 
in Pop B and RL AGN. 

VBLR Several lines of evidence (Corbin 1997a, b; Sulentic et al 2000b) suggest that a very broad 
line region (VBLR) of optically thin gas exposed to a very strong radiation field is located at the 
inner edge of the BLR. In the case of PG1416-129 the empirical evidence includes: (1) a difference 
of 9000 km/s between FWHM(BC) and FWHM (VBC) and (2) a much weaker response by the 
H/3vbc to a large continuum change. In an optically thick medium the intensity of a recombination 
line is governed by the luminosity of the ionizing continuum. If the medium is optically thin 
the intensity of the same recombination line is governed by the volume and density of the cloud 
distribution and is not directly related to the luminosity of the ionizing continuum. The much 
larger decline of H/3bc with respect to H/3vbc can therefore be explained if H/3bc is emitted in 
an optically thick medium, while a significant fraction ( ^ 50%) of H/3vbc is emitted by optically 
thin gas (see e.g. Shields et al. 1995). 

In the case of PG 1416-129, photoionization calculations performed using CLOUDY 94 (Fer- 
land 2000) are able to reproduce the H/3vbc luminosity assuming that an optically thin screen of 
gas surrounds the continuum sources at a distance of 10 17 ' 5 cm, with n e ~ 10 11 cm -3 , and T ^ 
10~ 2 . The H/3vbc is strongly redshifted and redward asymmetric (in PG 1416-129 the peak shift 
Av r ~ 850km s _1 , while the shift at line base is Av r ~ 4500 km s _1 ). If we ascribe the shift to 
gravitational plus transverse redshift (Corbin 1997b), then it is Av r ~ lSOMg^rj'g km s , and a 
black hole mass of a few 10 9 M is needed to produce Av r ~ 1000 km s~ 1 at r~10 17 ' 5 cm. 

Any optically thin screen with T ^ 0.01 is expected to be a strong source of HeiiA4686; for 
r ~ 0.1 HeiiA4686 can become comparable to H/3. In the case of PG 1416-129 milder conditions 
seem to be appropriate. We measure an upper limit of I(HeiiA4686vbc)/I(H/?vbc) ~ 0-25 (Table 
1), a value close to what has been held for a long time to be the canonical value for type 1 AGN. 
This is consistent with T ~ 10~ 2 deduced for the optically thin shell covering all the source needed 



to explain the luminosity of H/3vbc- 



3.3. Implication 3: An RL — RQ Population B Connection? 

We have recently described an Eigenvector 1 parameter space as the optimal discriminator 
between various AGN broad line classes (Sulentic et al. 2000a, b). Sources with FWHM (H/3bc) 
^ 4000 km s _1 (Population A) tend to be RQ while broader lined sources are much more often RL. 
We identified a RQ Population B that shows line profile and soft X-ray properties indistinguishable 
from, especially flat spectrum, RL sources. Most of the RL and 26% of the RQ sources in the PG 
sample (BG92) fall in the population B domain. The various planes of the correlation space show: 
(1) reasonably strong correlations among pop A sources and (2) little or no correlation among the 
RL and RQ POP B sources, although their mean values are an extension of the pop A correlations 
(Sulentic et al. 2000a). 

Observational commonalities (and at the same time differences from Pop A) between RL and 
RQ Pop B sources include: (i) stronger and more frequent optical variability (Ulrich et al. 1997), 
(ii) more complex/boxy/asymmetric profiles (Sulentic 1989; Eracleous & Halpern 1994; Marziani 
et al. 1996), (iii) occurrence of double-peaked profiles (e.g. Chen et al. 1989; Sulentic et al. 1995), 
(iv) occurrence of large single-peaked red/blue line shifts ( ^ 10 3 km s _1 ) and asymmetries (e.g. 
Marziani et al. 1993; Gaskell 1983), (v) the absence of a systematic ClvA1549 blueshift (Marziani 
et al. 1996; Sulentic et al. 2000a), (vi) the absence of a soft X-ray excess (Yuan et al. 1998; Sulentic 
2000a) and (vii) the presence of a VBLR emitting component in an uncertain number of sources. 
All line related comments except (v) refer to the Balmer lines where the phenomenology is better 
established. 

The VBLR property appears to be an important commonality between RL and RQ pop 
B. Sources like B2 1721-34/PKS1101-32 (Corbin 1997ab), PKS0837-12 (Corbin & Smith 2000), 
PKS0454-22 (Corbin & Boroson 1996), PKS0214+10 (Eracleous & Halpern 1994), PKS 0403-132/ 
PKS0405-123 (Marziani et al. 1996), 0159-117 (Brotherton 1996) are virtual twins of the 1990 
spectrum of PG 1416-129. Detection and accurate measurement of the VBLR feature requires at 
least moderate resolution ( ;$, 5A) and s/n ( ^ 20 in the continuum near H/3. RQ pop B analogs of 
PG1416-129 can be found in Marziani et al. (1996) where even relatively FEii opt strong RQ source 
Fairall 9 shows a VBLR feature. PG1416-129 shows us that pop B RQ sources can have a VBLR 
as strong as the strongest examples among RL AGN. 

The situation is less clear for much FEii op t stronger Population A sources. H/3bc is i n several 
cases observed to be blue-ward asymmetric, where the stronger blue wing may be associated with 
a high ionization wind emitting most of CivA1549 (Marziani et al. 1996; Sulentic et al. 2000ab). 
The available evidence suggests that RQ Pop A sources lack a VBLR like the one observed in Pop 
B objects (i.e. redshift~10 3 km/s and FWHM~10 4 km/s). Careful analysis of the Fen A4570- 
contaminated HeiiA4686 line profile would help to clarify the issue. 



3.4. Implication 4: Do We Always Measure the Same "BLR" in AGN? 

PG1416-129 suggests that the answer to this question is "no". Sources with "naked" VBLR 
lines certainly exist. A spectrum of PG1416-129 with lower s/n or resolution would not reveal the 
small BLR residual as a feature. In 1990 BG92 reported FWHM(H/3 B c) ~ 6100 km s" 1 and R Fen 
~ 0.2 for PG1416. The classical FWHM(H/3bc) was overestimated and RFeii underestimated. This 
mis-estimation of BLR properties will occur whenever a significant VBC is present unless the local 
continuum is set above that component. Using the 2000 spectrum to model the VBLR component in 
the 1990 spectrum (see Fig. 2 and Table 2) suggests that the correct values for the FWHM(H/?bc) 
and RfcII are 4000 km/s and 0.7 respectively. These parameter changes certainly exceed published 
error estimates for such measures. In the Eigenvector 1 correlation plane of FWHM(H/?bc) vs - 
RfcII, this change will move PG 1416-129 towards the population A domain. It will also produce 
an apparently FEll opt "stronger" source (from Rpcii~0.2 to 0.7) than is typically found for broad 
line RL/RQ pop B sources. 

If we parameterize PG1416-129 today we would measure FWHM ("H/3 S c") ~ 13000 km s" 1 
along with very weak/undetectable FEii opt . We expect no correlation between those two parameters 
because we know that FEII opt strength correlates with the classical H/3bc and not with the H/3vbc- 
We would in fact be measuring a "naked" VBLR. The implication for statistical studies and sample 
comparisons is clear: we must take into account the VBLR feature or we are not measuring the 
same thing in different sources. Sometimes we measure a pure BLR component, sometimes a pure 
VBLR component and sometimes a composite. The implications for the Eigenvector 1 space are 
also clear. It is possible that the pop A-B differences motivated by the PG sample reflect the 
presence of the extra VBLR emitting component in Pop B sources, and that the classical BLR may 
be more similar from source to source than we have appreciated. 



3.5. Implication 5: Inferences from Variability/Reverberation Studies 

Much of the predictive power of variability/reverberation studies involves determining the 
sequential response of the line core and wings. These results have previously been argued to 
support dominance of radial or rotational motions in the BLR (e.g. Gaskell 1988; O'Brien et al. 
1998; Goad et al. 1999). The implication of the VBLR is that the red wing of the line may have 
nothing to do with the line core. Therefore the presence or lack of a pre-, same or post response 
allows us to infer nothing about the source structure. Monitoring data for sources like Fairall 9 
could now be profitably reprocessed in the light of the VBLR concept. 

In a sample of 14 RL sources involved in an annual spectroscopic campaign Corbin & Smith 
(2000) report that most changes occur in the core of the line rather than the wings. They find that 
B2 1721-34 shows a decrease in the core component as the continuum declined. They also argue 
that the wings of the line have increased at the same time. Close examination of their spectra show 
that most of the change in the wings has occurred on the blue side. All of this is consistent with 



the idea that the peak is dominated by the classical BLR which is optically thick. The red wing, 
especially, arises in large part in optically thin gas, so it would be expected to change much less. 
In any case the two wings may be phenomenologically distinct. 

It is important to stress that AGN Pop A sources mirror the behavior of Pop B and RL AGN. 
The blue wing of H/3 seems to be associated with an optically thin wind emitting mostly HIL 
like CivA1549. The blue wing in population A sources may therefore show a lack of response to 
continuum change but for reasons unrelated to a VBLR component. 



4. Conclusion 

The quasar PG1416-129 shows two broad line regions. The classical BLR and VBLR are 
identified through inflections in the low ionization Balmer line profiles. They have a different place 
of origin and also arise from a physically distinct (opticaly thick/thin) regions as indicated by the 
difference of their response to the continuum change over the past ten years. 

MC, PM and JS acknowledge support from the Italian Ministry of University and Scientific 
and Technological Research (MURST) through grant Cofin 98-02-32. JS and TZ acnowledges 
support and telescope time from IA/UNAM. TZ acknowledges support from the Slovene Ministry 
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Fig. 1.— Spectra of the B./3 spectral region of PG 1416-129 on Feb. 17, 1990 and on May 8, 2000. 
Horizontal scale is rest-frame wavelength in A, vertical scale is observed specific flux (as in BG92). 
The two spectra have been normalized to the [Oiii]AA4959,5007 flux measured in 1990. 
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Fig. 2. — Continuum subtracted H/3 line profiles on 1990 Feb. 17 (upper panel) and on 2000 May 
8 (lower panel). The dotted line traces the original spectrum after FEll opt subtraction (FEii opt 
subtraction is however modest and has a negligible effect on the H/3 profile). The solid thin line 
represents the broad line spectrum, after narrow line and HeiiA4686 subtraction. The smooth, 
thick solid line in the lower panel is the result of a low order spline fitting to the Very Broad 
Component of H/3; in the upper panel, the fit result on the 2000 May 8 spectrum has been scaled 
by a factor 2.1 to match the 1990 Feb 17 observation. The H/3vbc spline fit has been scaled and 
shifted to account also for the HeiiA4686 profile (smooth thin lines). 
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Table 1. PG 1416-129 Emission Line Spectrum 3 



Flux b W c FWHM d 



[A] [k 



m s 



H7NC 


l.l c 


10: c 


< 400° 


H7bc+vbc 


11.4: f 


80: f 


> 4500 


[0 III]A4363 


2.2: c 


20: c 


1000: e 


He II A4686 N c 


0.6:§ 


1.0: s 


850 s 


He II A4686 V bc 


10.0 h 


75 h 


14000 h 


H/?NC 


2.0 


14 


250 


H/?BC+VBC 


40.0 


300 


9000 


[0 m]A4959 


9.0 


68 


370 


[0 III]A5007 


27.7 


205 


370 



a Line parameters measured on the 2000 May 8 spectrum unless otherwise noted, after rescaling 
to the [O III]A4959,5007 fluxes of the 1990 spectrum. 

Observed, in units of 10~" 15 ergs s _1 cm~ 2 . 

c Rest frame equivalent widths are measured against the bare continuum. 

Corrected for instrumental profiles assuming FHWMj nS ( r ~ 7 A and 3.5 A for the 1990 and 
2000 spectrum respectively. 

e Measured after subtraction of a scaled and shifted H/3bc+vbc to mimick H7BC+VBO 

f Measured scaling and shifting H/?bc+vbc- 
g Measured on the 1990 Feb. 17 spectrum. 

Measured using H/3vbc shifted and scaled to 1/4 its original flux. 

":" : uncertainty higher than standard estimate (see text). 
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Table 2. PG 1416-129 H/3 Emission Line Variability 



Line Identification 


1990 Feb. 17 a 




2000 May 8 b 




Flux 


W 


FWHM 


Flux 


W 


FWHM 






[A] 


[km s _1 ] 




[A] 


[km s _1 ] 


H/3bc+vbc 


86.0 


160 


6000 


40.0 


300 


9000 


H/3bc 


23.0 


47 


4000 


2.0 


13 


1450 


H/3vbc 


63.0 


110 


13000 


38.0 


220 


13000 



Specific flux at 4500 A F A « 7.5xl0~ 16 ergs s" 1 cm" 2 A -1 . 
b Specific flux at 4500 A F A ss 1.9X10" 16 ergs s" 1 cm" 2 A" 1 . 



